Disentangling perturbative and power corrections in precision tau decay analysis 
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Hadronic tau decay precision data are analyzed with account of both perturbative and power 
corrections of high orders within QCD. It is found that contributions of high order power corrections 
are essential for extracting a numerical value for the strange quark mass from the data on Cabibbo 
suppressed tau decays. We show that with inclusion of new five-loop perturbative corrections in the 
analysis the convergence of perturbation theory remains acceptable only for few low order moments. 
We obtain m s (M T ) — 130 ± 27 MeV in agreement with previous estimates. 

PACS numbers: 12.15.Lk, 13.35. Bv, 14.60.Ef 



Effects of strong interactions is a real stumbling block 
for investigating the electroweak sector of the Standard 
Model Q, 0, 0, EH • While there remain still many princi- 
pal problems of QCD as an underlying theory of strong 
interactions unresolved, an account of hadronic effects 
at the level of few percents is becoming a must for the 
high precision tests of the Standard Model and search 
for new physics d, El El El • Although the phenomenon 
of confinement is still beyond a complete quantitative 
theoretical explanation there is a solid qualitative under- 
standing of many features of QCD beyond perturbation 
theory that allows for a reliable use of perturbation the- 
ory (pQCD) in its applicability area for obtaining high 
precision predictions. The nonperturbative effects are 
accounted for through several phenomenological param- 
eters [Io|. A high precision achieved for the hadronic 
r-lepton decays both theoretically and experimentally 
makes the r-system a uniq ue testing ground of parti- 
cle interactions 0, 0, Ib'Sl ). The analysis of r-decays 
provides information usable in a variety of ways for: (i) 
extracting QCD parameters with high precision - strong 
coupling constant, s-quark mass, vacuum condensates of 
local operators within the operator product expansion; 
(ii) understanding general properties of perturbation the- 
ory and its asymptotic behavior at high orders; (iii) eval- 
uating the hadronic contributions necessary in the high 
precision tests of the Standard Model, e. g. the electro- 
magnetic coupling cxem(Mz), muon g — 2, Higgs mass. 

In this note we present a new analysis of the hadronic 
r-lepton decays with the main emphasize on the precision 
and reliability of the theoretical description within QCD. 

The total r-lepton decay rate into tau neutrino and 
hadrons normalized to the corresponding pure leptonic 
decay R T = T(t — > hv)/T(T — > Ivv) splits into a sum 
of strange and non-strange channels R T = i?f =0 + i?f =1 
with the experimental values Rr = 3.642 ± 0.012 and 
= 0.1625±0.0066 QEEIlB- The numerical values 
for the decay rates Rf =0 and i?f =1 are plausibly under- 
stood since in the parton model approximation the non- 
strange and strange parts of the decay rate are associated 
with the ud and us decay channels R^ =0 tx N c \V u d\ 2 



to the number of quark colors in QCD with 7V C = 3 
while the relative difference between the strange and 
non-strange channels is due to a numerical smallness of 
the V us entry of CKM matrix, \V US \ = 0.2196 ± 0.0026, 
as compared to the Cabibbo favored ud channel with 
\V ud \ = 0.9734 ± 0.0008 17]. For the determination of 
detailed characteristics of the spectrum in r-decays and 
further improvements upon precision the moments of the 
differential decay rate of the r lepton into hadrons 



R k J = 



have been extracted from the experimental data. The- 
oretically, the moments of the differential decay rate 
R kl are calculable in QCD perturbation theory (pQCD) 
within the operator product expansion for current corre- 
lators (OPE) 0. The range of indices (k,l) for the 
moments R kl in eq. Q should be properly chosen in 
order to guarantee the applicability of QCD perturba- 
tion theory for evaluating the moments with a strict 
control over the theoretical precision [l^. The analy- 
sis of the non-strange part of the decay rate results in 
an accurate determination of the strong interaction cou- 
pling constant a. A it) directly in the low energy domain 
for /i ~ M T |2oL I2H . The determination of the numer- 
ical value for the s-quark mass exploits the difference 
SR k T l = R kl s=0 /\V ud \ 2 - R kl s=1 /\V US \ 2 . This is a sensible 
setup as CKM is an external quantity to QCD and should 
be factor out. The analysis of the difference 6R kl is more 
demanding theoretically and requires much care in the 
interpretation of perturbation theory calculations in or- 
der to retain the full control over the obtained precision. 

The theoretical expression for the difference 5R kl is 
usually written in the general form 



6R k T l =2N c S EW J2 S n( m s,» s ) 

n>2 



(2) 



where Sew = 1-02 is the electroweak correction [22j 



The quantities 5„ 



,a s ) with n > 2 give corrections 



emerging within OPE technique in QCD. They include 
mass corrections within pQCD for n = 2 and power cor- 
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TABLE I: Perturbation theory moments of ud part 



(M) 


(Sp)1? 


(X \NLO 


/r \NNLO 
( d P)kl 


(0,0) 
(1,0) 
(2,0) 
(3, 0) 
(4,0) 


0.200 ± 0.005 
0.159 ±0.004 
0.135 ±0.004 
0.120 ±0.004 
0.110 ±0.004 


0.200 ±0.005 
0.167 ±0.006 
0.145 ±0.006 
0.137 ±0.007 
0.125 ±0.008 


0.200 ± 0.005 
0.167 ±0.007 
0.151 ±0.008 
0.144 ±0.009 
0.126 ±0.010 



5^(m s ,a s ) is very clean as they are computed through 
the two-point correlator of the weak charged currents 

i J dxe^(T3^x)3t(0))=q^n q (q 2 )+ giw n g (q 2 ) (3) 

with two scalar form factors II 9 (g 2 ) and H g (q 2 ). Such a 
decomposition of the general tensor correlator in eq. 
into a sum of scalar form factors avoids kinematical singu- 
larities. The function Tl g (q 2 ) receives contributions from 
the states with the total angular momentum J = 1 only 
that gives an attractive opportunity to analyze the decay 
data with respect to the spin content of the particles in 
the final state [2j|. The differential r-lepton decay rate 
is proportional to discontinuities i? g g (s) of the functions 
n 9iff (g 2 ) across the physical cut along the positive semi- 
axis q 2 — s > in the complex q 2 plane 



dR T 
ds 



oc 



M 2 ) 



(4) 



The numerical values for the moments of the ud part 
of the decay rate calculated in QCD are given in Ta- 
ble [I] The experimental input for the calculation is a 
perturbative correction dp to the total decay rate for the 
non-strange decays defined through 



N c Sew I Kd| 2 (l + 5p + Snp + Sew) 



with Sp = 0.200 ± 0.005. The power correction con- 
tribution Snp — —0.003 ± 0.004 is small and consis- 
tent with zero. The additive electroweak correction is 
also negligibly small, Sew — 0.001. This fixes the 
numerical value for the moment (A;, I) = (0, 0) which 
is used as the input for the determination of the nu- 
merical value for the strong coupling constant a s (M T ) 
from T-decays. The perturbation theory results are ob- 
tained in the approximation of massless u, d quarks and 
are rather stable up to the next-to-next-to-leading order 
(NNLO), i.e. including two corrections to the leading 
order non- vanishing result. The inclusion of the fifth or- 
der contribution in the coupling constant (NNNLO or 
three perturbation theory corrections to the leading or- 
der non-vanishing result with the "estimated" numerical 
value for the last coefficient of perturbation theory ex- 



still reasonable, we find (<5 P )f NNLO = 0.170 ± 0.007 and 
(<5p)30 NNL ° = 0.147±0.013. Thus, the pattern of conver- 
gence for the perturbation theory correction to the decay 
rate is 

(Mfo NNL ° = 0.170 = 0.159 ± 0.008 + 0.000 ± 0.003 
for the (1,0) moment and 



(Sp) 



NNNLO 
30 



= 0.147 = 0.120 ± 0.017 + 0.007 + 0.003 



for the (3, 0) moment. It is difficult to estimate the actual 
accuracy of the truncation of the asymptotic series in the 
coupling constant. 

The leading power corrections within the operator 
product expansion for the correlator from eq. are 
given by m 2 and quark condensate m s (ss) 



n«V)-n«V) = 



n"V)-n"V) 



47T 2 q 2 



.(5) 



The correction to H g (q 2 ) is ultraviolet divergent with /i 
being a subtraction point. 

The m 2 contribution to 5Rz l reads 



6?(a s ) = 3F klJ ^, 



Ft,° = S, 



/o 



(fc±2)m 



(0) 



The mass correction in the function H q (q 2 ) from eq. (0 
gives no contribution to the moments with I > in the 
leading order of perturbation theory |25j. The leading 
order result for the m 2 contribution from eq. 10 gets 
strongly renormalized in higher orders of perturbation 
theory [2(1 E7I Eij ] . In the approach based on the finite 
order perturbation theory one finds 



"'00 



5.333a s ± 46.0a 2 . ± (283.6 + k%)a 3 s ) (7) 



with a s — a s (M T )/ir. At the order there is an un- 
known constant fcf in the correlator II 9 (g 2 ) that con- 
tributes to the moments with I = 0. The use of the 
moments with I > allows for pushing the theoretical 
accuracy to the aj? level but the experimental data for 
such moments is less precise. For the numerical value of 
the strong coupling constant a s (M T ) = 0.344 ± 0.006 as 
extracted from the r lepton decay rate into non-strange 
hadrons the convergence of perturbation theory series is 
slow. The explicit convergence of higher moments within 
perturbation theory is worse, e.g. 
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^20 = •= (1 

5 v 



6.456a s ± 62.25a 2 + (547.8 + fcf )a,) . (8) 
The detailed analysis of convergence in the finite order 
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TABLE II: Coefficients of mi term 



TABLE III: Experimental moments 



(k,l) 


^kl 


pNLO 
?kl 


pNNLO 


(0, 0) 


1.877 ± 0.009 


2.913 ±0.040 


3.475 ± 0.067 


(1,0) 


1.997 ±0.016 


3.411 ±0.071 


4.408 ±0.134 


(2,0) 


2.135 ±0.023 


4.002 ±0.112 


5.613 ±0.237 


(3, 0) 


2.285 ±0.031 


4.689 ±0.165 


7.152 ±0.388 


(4,0) 


2.442 ± 0.040 


5.480 ±0.232 


9.109 ±0.609 



refs. 0, |2^. A new momentum to the implementa- 
tion of perturbation theory results to the analysis of r 
physics was given by developing the contour improved 
method for computation of the moments that allows to 
resum the effects of running of the cou pling co nstant 
in all orders of perturbation theory |29j, l30i I.'S1| . The 
technique is especially transparent in reformulation for 
the effective scheme description of the moments [22, |23| . 
Still the explicit convergence of the perturbation the- 
ory results for the m 2 correction is slower than for the 
ud part as is seen from Table ^ The next order gives 

F NNNLO = 4 on ± o 1Q2 and F NNNLO = U1Q8 ± 1 2gl 

with k\ — 160 |2s| . Note that recently the quantity 
analogous to k% has been computed for the diagonal vec- 
tor correlator 34] . The calculational techniques devel- 
oped for this evaluation (for details and further references 
see [34|) allows for computing the quantity k\ as well. 
The contribution of quark condensate 



rfcO 

o 4 



O) 



is linear in m s that implies a potentially large numeri- 
cal magnitude. It suffices to use the leading order ap- 
proximation of the coefficient function as the operator 
m s ss is renormalization group invariant that makes the 
change of the coefficient function small with running. In 
the numerical analysis we use (ss) — (0.8 ± 0.2) {uu) and 
(uu) = -(0.23 GeV) 3 H^. 
Further power corrections are written as 



n„ 



£ 

n>3 



{-q 2 ) n 



n 9 = E 

n>3 



2)(n-l) 



that leads to n > 3 contributions 

2 a 2n (fc + 2)! 



M T 2n (n- l)!(Jfe-n + 3)! 



(10) 



(11) 



with a 2n = {O q 2n ) - 2{0 9 2n+2 )/M 2 (see eq. @). Thus, 
the LO expression for 5R®° reads 



1 



■SR 



no 



ST 2 W !^ S) l7T 2 °* 



(12) 



2N C "" T M 2 M$ M«' 

The quantity ae = (Og) — 2 (Of) /M 2 contains a contribu- 



(k,l) 


(0,0) 


(1,0) 


(2,0) 


(ffl^) exp 


0.394 ±0.137 


0.383 ± 0.078 


0.373 ± 0.054 



g-part of the correlator @. Note that this part receives 
only J — 1 contributions. 

The perturbation theory coefficients Fu of the m 2 con- 
tribution are given in Table UTI up to NNLO. An account 
of new NNNLO results from ref. [24\ shows that the con- 
vergence virtually disappears for k > 2. Thus, the (2,0) 
moment is indeed marginal for the perturbation theory 
calculations to be trusted. The quantity SR^ a is sensitive 
to the power corrections up to a 2 k+G which include the 
maximal dimension operator (O^+g)- 

Experimental results for the moments of the decay rate 
from ref. 0] are given in Table ITTT1 Recently published 
results are rather close 36] . Using these data and the- 
oretical calculation at the next-to-next to leading order 
(NNLO) of perturbation theory we extract the numerical 
value for m s . With ae — 0.001 GeV 6 which is obtained 
in the approximation (Of) = and neglecting the power 
corrections a 2n with n > 3 one finds the following results 
at NNLO: 



(0,0) 
(1,0) 
(2,0) 



m s (M 2 ) 
m s {M 2 ) 
m s {M 2 ) 



130 ± 27 oxp ± x t h MeV 
110±13 exp ±:rth MeV (13) 
94 ± 8 oxp ± xth MeV 



Here xth denotes a theoretical uncertainty of a given mo- 
ment to be discussed in much detail later. Note that the 
error in m s due to the experimental uncertainty of the 
measured moments is large. Thus the linear approxima- 
tion (propagation of errors) is not really applicable for the 
error analysis of the extracted value of m s . Considering 
the leading order of perturbation theory for the illus- 
tration we obtain the following values for the extracted 
masses 



(0,0) 
(1,0) 
(2,0) 



m s {M 2 )\ LO 
m s {M 2 )\ LO 
m s (M 2 )\ 



LO = 



171 ± 37 exp MeV 

155 ± 19 oxp MeV (14) 

141 ± 13ox P MeV. 



For the low value 0.394 - 0.137 = 0.257 of the moment 
(0, 0) we find m s {M 2 ) = 171 - 37 = 134 MeV while for 
the high value 0.394 ± 0.137 = 0.531 we have m s {M 2 ) 
1 7 1 ± 31 = 202 MeV. In our analysis we take the biggest 
error as a conservative estimate for the uncertainty of the 
extracted numerical value for the quark mass. 

It is instructive to analyze how the theoretical predic- 
tions for the moments are composed from the different 
QCD contributions. To see the relation between pertur- 
bation theory and nonperturbative contributions is most 
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the numerical values for the experimental moments into 
the contributions of perturbation theory and power cor- 
rections at NNLO: 



(0,0) 
(1,0) 
(2,0) 



0.394 = 0.334 + 0.060 
0.383 = 0.306 + 0.077 
0.373 = 0.286 + 0.087 



(15) 



The first term is a "trivial" perturbation theory power 
correction term proportional to the mass squared of the 



strange quark 



-term. It is explicitly independent of 



the vacuum structure of QCD or the numerical values 
of the vacuum condensates. The second term is given 
by the quark condensate and is linear in m s that makes 
it potentially large numerically. This term is renormal- 
ization group invariant. Note also that by construction 
the ratios of quark condensate contributions in different 
perturbation theory orders are equal to the ratios of the 
corresponding values of quark masses. One sees that the 
contribution of the leading power correction is numeri- 
cally significant. The relative magnitude of the quark 
condensate contribution increases for larger (k, 0) mo- 
ments from 18% to 30%. This calls for the analysis of 
contributions of higher order power corrections. 

There are no general systematic techniques to estimate 
the theoretical errors due to truncation of the asymptotic 
series. This is really difficult problem that should be 
considered for any given case. Let us consider a pattern 
of convergence for the extracted masses in different orders 
of perturbation theory for coefficient functions: 



(0,0) 
(1,0) 
(2,0) 



m s (M T 2 ) 
m s (M 2 ) 
m s (M T 2 ) 



130 = 171 - 30 - 11 MeV 
110 = 155 - 31 - 14 MeV (16) 
94 = 141 - 32 - 15 MeV. 



Taking one half of the last term as the estimate of the 
truncation error we find the following results 



(0,0) 
(1,0) 
(2,0) 



m s (M 2 ) = 130 + 6 MeV 
m s (M 2 ) = 110 + 7 MeV 



(17) 



m s (M 2 ) 



94 + 8 MeV. 



This uncertainty comes from the truncation of the pertur- 
bation theory series for the coefficients F^i- Combining 
them with the experimental errors we finally get 



(0,0) 
(1,0) 
(2,0) 



m s (Af 2 ) = 130 ± 27 oxp ± 6 tr MeV 
m s (M 2 ) = 110 ± 13 cxp ± 7 tr MeV (18) 
m s (Af 2 ) 



cxp 

94 ± 8 cxp ± 8 tr MeV. 



For higher moments the relative weight of truncation er- 
ror is larger. 

Still, this is not a whole story and we have to esti- 
mate uncertainty due to higher power corrections. One 
possibility to perform such an estimate is to consider the 



there is no contribution of higher power correction terms 
that leaves it untouched. 

For the (1,0) moment there is a contribution of one 
higher power correction. Taking again one half of the 
last term we get the uncertainty 0.077/2 w 0.04 from 
eq. I|15fl ■ This uncertainty leads to an additional error in 
the mass of 7 MeV (this number can be obtained either 
by a direct computation or from the results for the ex- 
perimental error in the linear approximation, it should 
be around one half of it, 13/2 k, 7). 

For the (2, 0) moment there are at least two next power 
correction contributions. This fact can make the result- 
ing error larger if these corrections are correlated. Still 
taking one half of the last term we get the uncertainty 
0.087/2 = 0.044 from eq. EH This uncertainty leads to 
the additional error in the extracted mass of 7 MeV. Hav- 
ing in mind the possibility of correlation of power correc- 
tions one could add 50% to the error (see also eq. (|20l) ) 
and finally get 10 MeV. 

This analysis leads to the following uncertainties due 
to power corrections that should be added to the pertur- 
bation theory uncertainties from eq. I|17(l 



(0,0) 
(1,0) 
(2,0) 



m s (A/ 2 ) = 130 + 6 + MeV 
m s (A/ 2 ) = 110 + 7 + 7 MeV (19) 
m s (A/ 2 ) = 94 + 8 + 10 MeV. 



These uncertainties are only indicative as it is really dif- 
ficult to make any reliable evaluation of the errors due to 
higher power corrections. 

A remarkable feature of the obtained results for the nu- 
merical value of the strange quark mass is the strong de- 
pendence of the extracted numerical value for m . {Mz) on 
the particular moment used for the analysis 32, 371 l38j . 
The theoretical errors are dominated by the uncertainty 
due to the truncation of the perturbation theory series. 
The experimental and theoretical errors are independent 
and added in quadrature give the total uncertainty of 
28, 16, 15 MeV for the m s values extracted from the 
(0,0), (1,0), (2,0) moments. From perturbation theory 
point of view the moment (0, 0) is the best one but its 
experimental accuracy is low. With neglecting the contri- 
butions due to higher power corrections the value of m s 
extracted from the (1,0) moment has a smaller total un- 
certainty because the experimental accuracy for the (1,0) 
moment is much better than that of the (0, 0) moment. 
This analysis makes clear that in order to improve upon 
the total precision one should reduce the experimental 
errors for the low moments. 

Three estimates for the numerical value of m s in 
eq. i|18|) are consistent within the error bars since the 
experimental errors are rather large. However, the sit- 
uation is definitely unsatisfactory as there is a system- 
atic decrease of the central value for m s extracted from 
different moments. The reason for this systematic de- 
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tions in the analysis. The analysis of experimental data 
with the inclusion of higher power corrections as free pa- 
rameters of the fit gives m s (M 2 ) = 130±27 eX p±6th MeV, 
a 8 = 0.05 GeV 8 , am = -0.3 GeV 10 . The value for m s 
has not changed since the higher power corrections do 
not contribute to the (0, 0) moment. However the equa- 
tions for moments (0, 1) and (0, 2) are satisfied with the 
same m s because of contributions due to higher power 
terms as.io- It is instructive to see how the moments are 
composed with an account of the contributions of power 
corrections (ordered by the dimensionality) 

(0,0): 0.394 = 0.341| m 2 + 0.061| ms - 0.008| ae (20) 
(1,0): 0.383 = 0.433| m | +0.092| ms -0.024| O6 
-0.118| O8 

(2,0): 0.373 = 0.552 | m j + 0.123 | TOs - 0.048 1„ 6 
-0.472| Q8 +0.218| aio . 

There is a sizable contribution from higher power correc- 
tions to the high moments. The numerical values for the 
high dimensional condensates are extracted as fit param- 
eters from the data. Can they be understood from the 
present knowledge of the numerical magnitude of power 
corrections? The structure of power corrections for the 
relevant correlator is a^n = C2„m s A 2 "' 1 with A being 
a nonperturbative infrared scale of QCD and the fac- 
tor m s appears as the power corrections to the differ- 
ence of the correlators should vanish in the chiral limit. 
For m s = 130 MeV and A ~ m p = 0.77 GeV one finds 
a 8 = 0.02Cg and am = 0.012Cio. One expects a fast 
growth of coefficients C2n with n, e.g. |25j. Assuming 
a factorial growth Cm ~ fi\ one obtains for the ratio 
aw/ as = 5A 2 = 3 GeV 2 while the values extracted from 
the data is 6 GeV 2 . This is quite reasonable result given 
the simplicity of the estimate. Note that the value of 
ae seems to be very small within this logic. It is known 
however that there is a huge numerical cancellation be- 
tween first several low dimensional power corrections in 
the vector and axial channels. Our method of estima- 
tion cannot catch such delicate features of the physical 
spectrum. 

One can also estimate the magnitude of power correc- 
tions phenomenologically. Assuming for simplicity that 
the continuum contribution basically cancels in the dif- 
ference of ud and us channel one is left with the difference 
of the resonance contributions. For the axial part it is 
given by the difference 



n«V)-iW 



f 2 

J ud 



f, 



2 

us 



(21) 



' us 



the corresponding channels. Taking f£ d 



f, 



2 

us 



0.6/4tt 2 GeV 2 and m ud 



m Kl 



1.4 GeV one finds |ag| 



= 1.24 GeV, m us 
0.06 GeV 8 and [oiol 



fit. The contribution of the vector channel to the differ- 
ence is much smaller because the masses of corresponding 
resonances p(770) and i^* (890) are smaller. 

It seems hopeless to estimate numerical values for high 
dimension contributions from first principles. Even if the 
basis of relevant operators is identified and the coefficient 
functions are found one remains with the problem of nu- 
merical values for vacuum condensates. The factorization 
hypothesis becomes less trustworthy with increasing di- 
mension of the condensates. Still our consideration shows 
that the fit to data is consistent with the estimates based 
on the structure of hadronic spectra. The emerging pic- 
ture is fairly sensible from the physical point of view. 
The large fc-moments give a high energy resolution of the 
spectrum as they are saturated with quite a limited part 
of the total spectrum; in fact, the large fc-moments virtu- 
ally become the exclusive observables. For such almost 
exclusive observables the influence of high order power 
corrections should be large. 

Note that the precision of experimental data is not 
sufficiently good that allows for a (marginal) fit of all the 
moments within the error bars with only one value of 
the strange quark mass and without including the higher 
power corrections. We do not discuss this possibility as 
power corrections are known to be definitely included as 
free parameters of the fit. 

The obtained numerical value for the strange quark 
mass is compatible with results available in the liter- 
ature 0, 0, EH E3- However the central value is 
larger than the lattice results 0, 0, 0|. We obtain 
m s (2 GeV) = 125 ± 28 MeV. 

To conclude, our current analysis of the r-decay data 
shows that the requirement of the applicability of per- 
turbation theory strictly limits the allowed range of mo- 
ments with the (k,l) — (2,0) moment being almost 
marginal since for the higher moments the convergence 
of perturbation theory series is virtually absent. The 
contributions of power corrections are large for the high 
fc-moments and have to be retained that stabilizes the 
value of the strange quark mass extracted from the data. 
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